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It is a challenge to manage complex systems efficiently without confronting NP-hard problems. 
To address the situation we suggest to use self-organization processes of prime integer relations 
for information processing. Self-organization processes of prime integer relations define correlation 
structures of a complex system and can be equivalently represented by transformations of two- 
dimensional geometrical patterns determining the dynamics of the system and revealing its structural 
complexity. Computational experiments raise the possibility of an optimality condition of complex 
systems presenting the structural complexity of a system as a key to its optimization. From this 
perspective the optimization of a system could be all about the control of the structural complexity 
of the system to make it consistent with the structural complexity of the problem. The experiments 
also indicate that the performance of a complex system may behave as a concave function of the 
structural complexity. Therefore, once the structural complexity could be controlled as a single 
entity, the optimization of a complex system would be potentially reduced to a one-dimensional 
concave optimization irrespective of the number of variables involved its description. This might 
open a way to a new type of information processing for efficient management of complex systems. 

PACS numbers: 89.75.-k, 89.75.Fb 



I. INTRODUCTION 

It is a challenge to manage complex systems efficiently 
without confronting NP-hard problems. To address the 
situation we consider the description of complex systems 
in terms of self-organization processes of prime integer 
relations [l[ and suggest to use the processes for infor- 
mation processing. 



II. THE HIERARCHICAL NETWORK OF 
PRIME INTEGER RELATIONS 

The description is realized through the unity of two 
equivalent forms, i.e., arithmetical and geometrical. We 
briefly present the forms in order to introduce the hier- 
archical network of prime integer relations. More details 
may be found in [l| , 0] . 

A. THE ARITHMETICAL FORM 

In the arithmetical form a complex system is char- 
acterized by hierarchical correlation structures built in 
accordance with self-organization processes of prime in- 
teger relations. As each of the correlation structures is 
ready to exercise its own scenario and there is no mecha- 
nism specifying which of them is going to take place, an 
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intrinsic uncertainty about the complex system exists. 
At the same time, the information about the correlation 
structures can be used to evaluate the probability of an 
observable to take each of the measurement outcomes. 
Therefore, the arithmetical form of the description pro- 
vides the statistical information about a complex system. 

The form reveals nonlocal correlations without refer- 
ence to signals as well as the distances and local times 
of the parts. Thus, the arithmetical form suggests that 
parts of a complex systems may be far apart in space and 
time and yet remain interconnected with instantaneous 
effect on each other, but no signalling. Namely, if a cor- 
relation structure of a system is selected and some parts 
are specified, then through the prime integer relations in 
control of the correlation structure the other parts are 
immediately defined. 

Through the arithmetical form of the description we 
become aware of the hierarchical network, i.e., a set 
of mutually self-consistent elements built by all self- 
organization processes of prime integer relations. Arith- 
metic ensures that not even a minor change can be made 
to any element of the network (Figure 1). It is helpful 
to think that self-organization processes of prime integer 
relations take place in the hierarchical network as they 
compose more and more complex structures. 

The hierarchical network of prime integer relations is 
a causal structure. As parts of a system change and the 
prime integer relations in control of the system cause the 
other parts to change accordingly, an event takes place. 
Once the changes have been realized, the event is fixed 
in space and time with respect to the reference frames 
of the parts. For the parts the effect of the event has 
not necessarily be the same, but for each part it is ap- 
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FIG. 1: The figure shows a hierarchical structure of prime 
integer relations defined as a result of a self-organization pro- 
cess starting at the zero level with integers 16, 13, 11, 10, 7, 6, 
4, 1 in the "positive state" and integers 15, 14, 12, 9, 8, 5, 3, 2 
in the "negative" . The process is controlled by arithmetic and 
cannot progress to level 5, because arithmetic determines that 
+16 4 -15 4 -14 4 + 13 4 -12 4 + ll 4 + 10 4 -9 4 -8 4 + 7 4 + 6 4 -5 4 + 
4 4 _ 3 4 _ 2 4 + 1 4 7^ 0. This hierarchical structure can specify a 
correlation structure of a system made of 16 elementary parts 
Pi,i — 1, 16 whose changes Axi, i = 1, 16 of the observ- 
ables Xi,i — 1, 16 relative to their reference frames can be 
given by the sequence Axi...Aa;i6 =+1 — 1 — 1 + 1 — 1 + 1 + 
1-1-1 + 1 + 1-1 + 1-1-1 + 1. 



propriately determined by the prime integer relations. 
However, the prime integer relations at work for a sys- 
tem have no causal power to effect systems controlled by 
separate prime integer relations. As a result, informa- 
tion about the systems is blocked for the observers of the 
system. 



B. THE GEOMETRICAL FORM 

Specified by two parameters e > and S > the ge- 
ometrical form arises as the self-organization processes 
of prime integer relations find isomorphic realization in 
terms of transformations of two-dimensional geometrical 
patterns Q. As a result, hierarchical structures of prime 
integer relations defining the correlation structures of a 
complex system become equivalently represented by hi- 
erarchical structures of geometrical patterns determining 
the dynamics of the system and revealing its complexity. 
The quantitative description of the system turns out to 
be given by the description of the geometrical patterns 



Figure 2 shows a hierarchical structure of geometrical 
patterns, which for given e > and 5 > is isomorphic 
to the hierarchical structure of prime integer relations 
depicted in Figure 1. A scale invariant property and a 
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FIG. 2: The figure shows a hierarchical structure of geomet- 
rical patterns, which for given e and 5 is isomorphic to the 
hierarchical structure of prime integer relations shown in Fig- 
ure 1. Under the integration of the function ^/j', *?i +1 '(to) = 
0, / = 0,1,2,3 the geometrical patterns of the level I form 
the geometrical patterns of the higher level I + 1. The width 
Wi of a geometrical pattern at level I, I = 1, ...,4 equals the 
sum Wi = Wi-i } ieft + Wi-i } right = 2Wi-i of the widths of 
the two geometrical patterns it is made up of, where Wi-\^ e f t 
and Wi~\ t ri g ht are the widths of the left and right geometrical 
patterns and Wi-i,i e ft = Wi-i, r i g ht = Wi-i as each geomet- 
rical pattern at level / — 1 has the same width Wi-.± = 2 l ~ 1 e 
and ti = ei,i = 0, 1, 16. The height Hi of the geometrical 
pattern equals Hi — where Si-i is the area of each of the 

geometrical patterns it is composed of. The width Wi and the 
height Hi can specify length scales of the level I, I = 0, 1, 4. 



renormalization group transformation come to our atten- 
tion, as we consider the connection between a geometrical 
pattern and the geometrical patterns it is made of. 

Although the geometrical patterns are not triangles at 
level I, I = 2, 3, 4, yet the boundary curve, as the graph of 
the function w£ , is such that the area Si of a geometrical 
pattern can be simply given, as if it were a triangle, by 
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where W\ and Hi are its width and height. 

The renormalization group transformation at level 4 
defines 
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and uses a coarse-grained procedure replacing the geo- 
metrical pattern made up of 8 geometrical patterns at 
level 1 by their enlarged version with the boundary curve 
as the graph of the function ■ Each of the geometrical 
patterns at level 1 is elementary in the sense that it is 
fully specified by a prime integer relation made up of two 
integers without further internal structure. The width 
2e' and the height S'e' of the renormalized geometrical 
pattern at level 4 are given in terms of e' and 6' in the 
same way as the width 2s and the height Se of a geomet- 
rical pattern at level 1 are given in terms of e and S. The 
two geometrical patterns at level 4 have the same width, 
height and area 

/ V [ * ] (t)dt = / "f l 2 ] {t)dt, 

J t J t 

but the lengths of their boundary curves are different 
(Figure 2). 

The arithmetical and geometrical forms unite the dy- 
namics and the structure of a complex system as two 
sides of the same entity in the preservation of the system 
as a whole. In particular, at one side the dynamics of 
the parts is determined to produce spacetime patterns of 
the parts to fit precisely into the geometrical patterns of 
the complex system. Under this condition at the side the 
corresponding prime integer relations can provide the re- 
lationships between the parts for the complex system to 
exist. If the spacetime patterns do not fit even slightly, 
then one or more of the relationships are not in place and 
the complex system collapses. 

To measure the complexity of a system in terms of self- 
organization processes of prime integer relations a con- 
cept of structural complexity is introduced Starting 
with the integers at the zero level, the self-organization 
processes of prime integer relations progress to different 
levels and thus produce a hierarchical complexity order. 
In particular, the higher the level self-organization pro- 
cesses progress to, the greater is the structural complex- 
ity of a corresponding system. In our description sys- 
tems can be compared in terms of complexity by using 
two equivalent forms: in structure - by the hierarchical 
structures of prime integer relations and in dynamics - 
by the hierarchical structures of geometrical patterns. 

Remarkably, based on integers and controlled by arith- 
metic only self-organization processes of prime integer re- 
lations can describe complex systems by information not 
requiring further simplification. 

III. TESTING A NEW MEDIUM FOR 
INFORMATION PROCESSING 

The correlation structures of a complex system contain 
information about the parts. By changing some parts the 
information can be processed as the other parts change in 
accordance with the correlation structures. This shows 
the importance of self-organization processes of prime in- 
teger relations for information processing. Namely, for a 



given problem they could be used to build the correlation 
structures of a system in processing information demon- 
strating the optimal performance for the problem. 

As a result, we suggest the hierarchical network of 
prime integer relations as a new medium for information 
processing and are interested in its navigating proper- 
ties. It would be important if for any problem the per- 
formance of a system could behave as a concave function 
of its structural complexity. Guided by this property the 
performance global maximum for a problem could be ef- 
ficiently found. It would be also beneficial if at the global 
maximum the structural complexities of the system and 
the problem could be related through an optimality con- 
dition. 

The optimality condition might be interpreted as fol- 
lows: if through arithmetical interdependencies emerging 
in the hierarchical network between a computing system 
and a problem a new building block is formed at the high- 
est possible level, then the optimal performance takes 
place. Or, in other terms, a computing system finds the 
solution to a problem, once arithmetic interdependencies 
emerging between the system and the problem provide a 
channel to obtain the desired information. 

It is worth to note that since the correlation structures 
of a system are completely determined by prime integer 
relations, which are equivalent to two-dimensional geo- 
metrical patterns, the entropy of the system, measuring 
its information content, can be connected with the areas 
of the two-dimensional patterns. Thus, in our approach 
there is a general connection between entropy and area. 

Computational experiments have been conducted Q 
to test the navigating properties. In particular, an opti- 
mization algorithm A, as a complex system, of N com- 
putational agents minimizing the average distance in the 
travelling salesman problem (TSP) is developed. The 
agents work in parallel and start in the same city by 
choosing the next city at random. Then an agent at each 
step visits the next city by using one of the two strategies: 
a random strategy or the greedy strategy. 

In the solution of a problem with n cities the state of 
the agents at step j,j — 1, n — 1 can be described by a 
binary sequence Sj — Sij...SNj, where = +1, if agent 
i,i = 1,...,N uses the random strategy and sy = — 1, 
if the agent i uses the greedy strategy. The dynamics 
of the complex system is realized as the agents step by 
step choose their strategies and can be encoded by an 
N x (n — 1) binary strategy matrix 

S = {sij,i = 1, ...,N,j = 1, ...,n - 1}. 

We try to change the structural complexity of the al- 
gorithm A monotonically by forcing the system to make 
the transition from regular behaviour to chaos by period- 
doubling. To control the system in this transition a pa- 
rameter v, < v < 1 is introduced. It specifies a thresh- 
old point dividing the interval of current distances passed 
by the agents into two parts, i.e., successful and unsuc- 
cessful. This information is provided for an optimal if- 
then rule |l[ that each agent uses to choose the next strat- 
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egy. The rule relies on the Prouhet-Thue-Morse (PTM) 
sequence 

+ 1-1-1 + 1-1 + 1 + 1-1 ... 

and has the following description: 

1. if the last strategy is successful, continue with the 
same strategy. 

2. if the last strategy is unsuccessful, consult PTM 
generator which strategy to use next. 

Remarkably, for each problem p tested from a class 
V it has been found that the performance of the algo- 
rithm A indeed behaves as a concave function of the 
control parameter with the only global maximum at a 
value v*(p). The global maximums {v*(p),p £ V} are of 
interest to probe whether the structural complexities of 
the algorithm A and the problem are related through an 
optimality condition. For this purpose strategy matrices 

{S(v*(p)), P eV} 

corresponding to the global maximums {v*(p),p G V} 
are used and the structural complexities of the algorithm 
A and a problem p are approximated as follows. The 
structural complexity C(A(p)) of the algorithm A is ap- 
proximated by the quadratic trace 

1 1 N 

C(A{p)) = w tr(V 2 (v*{p))) = W T, X i 

of the variance-covariance matrix V(v* (p) ) obtained from 
the strategy matrix S(v*(p)), where Xi,i — 1,...,N are 
the eigenvalues of V(v*(p)). 

The structural complexity C(p) of the problem p is 
approximated by the quadratic trace 

1 " 

C(p) = -tr(M*(p)) = -Y;(K? 

i—l 

of the normalized distance matrix 

M(p) = {dij/dmax^J = l,...,n}, 

where X^i — 1, ...,n are the eigenvalues of M(p), dij 
is the distance between cities i and j and d max is the 
maximum of the distances. 

To reveal the optimality condition the points with the 
coordinates 

{x = C(p),y = C(A(p)), P eP} 

are considered. The result has been indicative of a linear 
relationship between the structural complexities and thus 
suggests an optimality condition of the algorithm A Q : 

// the algorithm A demonstrates the optimal perfor- 
mance for a problem p, then the structural complexity 
C(A(p)) of the algorithm A is in the linear relationship 

C{A{p)) = 0.67C(p) + 0.33 



with the structural complexity C'(p) of the problem p. 

According to the optimality condition if the optimal 
performance takes place, then in terms of the structural 
complexity the dynamics of the algorithm A is in a cer- 
tain relation with the structure of the problem p, i.e., the 
distance network with the vertices as the cities and the 
edges specifying the pairwise distances. 

The optimality condition is a practical tool. For a 
given problem p by using the distance matrix we can 
calculate the structural complexity C(p) of the problem 
p and from the optimality condition find the structural 
complexity 

C(A{p)) = 0.67C(p) + 0.33. 

Then to obtain the optimal performance of the algorithm 
A for the problem p we need to tune the control param- 
eter for the algorithm A to function with the structural 
complexity C(A(p)). 



IV. THE ORDER OF THE PROCESSES AND 
EFFICIENT QUANTUM ALGORITHMS 

The computational results point that by using self- 
organization processes of prime integer relations it may 
be possible to design classical algorithms comparable to 
efficient quantum algorithms. 

Quantum algorithms rely on the practical use of entan- 
glement, whose sensitivity challenges the development of 
relevant technologies. In a TSP quantum algorithm the 
wave function would be evolved to maximize through the 
amplitudes the probability of the shortest routes to be 
measured. However, there is no general direction known 
for the evolution to take in order to make the quantum 
algorithm efficient. In this regard the majorization prin- 
ciple seems to play an important role [4J. It provides 
a local navigation, but without information about the 
whole landscape. 

While the nature of quantum entanglement is yet to be 
understood [j| in our approach the nonlocal correlations 
are known from their origin in the self-organization pro- 
cesses of prime integer relations. The question is whether 
this knowledge could be used to provide computational 
resources comparable to quantum computation. 

In this paper we focus on whether such a resource could 
be obtained from the nonlocal correlations although used 
in classical computation, but with the order of the pro- 
cesses preserved. If that could be possible, then the or- 
der of the self-organization processes would establish a 
general direction for efficient computation. Remarkably, 
the experiments raise the possibility that if the evolution 
goes in this direction, then the performance landscape 
becomes concave. 

In the experiments we use the parameter to control 
the correlation structures of the computing system with 
their consequences observed in the routes taken by the 
agents. To help in associations with the quantum case 
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the average distance for a value v of the parameter can 
be written as 

D = ^( 7 i,...,n-i(«M[l ) ... J n-l >) + ... 

+7n-i,...,i(v)d([n - >)) 

where 7»i,...,t„_i(«) is the number of the agents followed 
the route [ii, ...,i n -x >, i n -i >) is its distance 

and the n cities are labelled by 0, 1, ...,n — 1 with the 
initial city by 0. The interpretation of the coefficient 

N 

as the probability of the route [ii, > suggests that 
the minimization of the average distance considered in 
the algorithm A can be connected with the maximization 
of the probability of the shortest route considered in TSP 
quantum algorithms. 

Significantly, the experiments have shown that in the 
case of the algorithm A the maximization of the prob- 
ability turns out to be a one-dimensional concave opti- 
mization. In its course the computing system adapting 
to the problem got more complex or simpler until the 
global maximum is reached and the structural complex- 
ities of the computing system and the problem become 
related through the optimality condition. More connec- 
tions might arise if the wave function could be involved 
in the description of the correlation structures. 

We note that the algorithm A shares a common fea- 
ture with Shor's algorithm, which also relies on the PTM 
sequence @. 

V. CONCLUSIONS 

We have suggested that self-organization processes of 
prime integer relations could be used for information pro- 



cessing. 

In particular, for a given problem self-organization pro- 
cesses of prime integer relations could be used to build 
the correlation structures of a system in processing infor- 
mation demonstrating the optimal performance for the 
problem. Remarkably, the processes can be equivalently 
represented by transformations of two-dimensional geo- 
metrical patterns determining the dynamics of the sys- 
tem and revealing in the information processing its struc- 
tural complexity. 

The information processing would be distinctive, be- 
cause self-organization processes of prime integer rela- 
tions can define the correlation structures of a system 
without reference to the distances, local times and sig- 
nals between the parts. 

Computational experiments testing competitive ad- 
vantages of the information processing have been pre- 
sented. They raise the possibility of an optimality condi- 
tion of complex systems: if the structural complexity of 
a system is in a certain relationship with the structural 
complexity of a problem, then the system demonstrates 
the optimal performance for the problem. 

The optimality condition presents the structural com- 
plexity of a system as a key to its optimization. From 
its perspective the optimization of a system could be all 
about the control of the structural complexity of the sys- 
tem to make it consistent with the structural complexity 
of the problem. 

Importantly, the experiments also indicate that the 
performance of a complex system may behave as a con- 
cave function of the structural complexity. Therefore, 
once the structural complexity could be controlled as a 
single entity, the optimization of a complex system would 
be potentially reduced to a one-dimensional concave opti- 
mization irrespective of the number of variables involved 
its description. This might open a way to a new type of 
information processing for efficient management of com- 
plex systems. 
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